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A method of finding general solutions of second-order nonlinear ordinary differen- 
tial equations by extending the Prelle-Singer (PS) method is briefly discussed. We 
explore integrating factors, integrals of motion and the general solution associated 
with several dynamical systems discussed in the current literature by employing our 
modifications and extensions of the PS method. In addition to the above we intro- 
duce a novel way of deriving linearizing transformations from the first integrals to 
linearize the second order nonlinear ordinary differential equations to free particle 
equation. We illustrate the theory with several potentially important examples and 
show that our procedure is widely applicable. 
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1. Introduction 

Solving nonlinear ordinary differential equations (ODEs) is one of the classical 
but potentially important areas of research in the theory of dynamical systems 
(Arnold 1978; Jose & Saletan 2002; Wiggins 2003). Indeed the last century has wit- 
nessed considerable amount of research activity in this field. Progress has been made 
broadly in two different ways: one is through geometrical analysis and the other is 
through analytical studies. The modern geometrical theory was born with Poincare 
and developed later vigorously by Arnold, Moser, Birkhoff and others (Percival & 
Richards 1982; Guckenheimer & Holmes 1983; Wiggins 2003). In parallel with the 
geometrical theories various analytical methods have also been devised to tackle 
nonlinear ODEs. The ideas developed by Kovalevskaya, Painleve and his coworkers 
have been used to integrate a class of nonlinear ODEs and obtain their underly- 
ing solutions (Ince 1956). As a consequence of these studies, nonhnear dynamical 
systems are broadly classified into two categories, namely, (i) integrable and (ii) 
nonintegrable systems. Indeed one of the important current problems in nonlinear 
dynamics is to identify integrable dynamical systems (Lakshmanan & Rajasekar 
2003). Of course, these methods have close connection with the group theoretical 
approach introduced by Sophus Lie in the nineteenth century and subsequently ex- 
tended by Cartan and Tresse, to integrate ordinary and partial differential equations 
(see for example Olver (1995); Bluman & Anco (2002)). 

In order to identify such integrable dynamical systems different techniques have 
been proposed, including Painleve analysis (Conte 1999), Lie symmetry analysis 
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(Bluman & Anco 2002) and direct methods of finding involutive integrals of motion 
(Hietarinta 1987). Each method has its own advantages and disadvantages. For a 
detailed discussion about the underlying theory of each method, their limitations 
and applications we refer to Lakshmanan & Rajasekar (2003). Also certain nonlinear 
ODEs can be solved by transforming them to linear ODEs whose solutions are 
known. In fact, linearization of given nonlinear ODEs is one of the classical problems 
in the theory of ODEs whose origin dates back to Cartan. For the recent progress 
in this direction we refer Olver (1995). 

In this direction sometime ago Prellc & Singer (1983) have proposed a proce- 
dure for solving first order ODEs that presents the solution in terms of elementary 
functions if such a solution exists. The attractiveness of the PS method is that if 
the given system of first order ODEs has a solution in terms of elementary functions 
then the method guarantees that this solution will be found. Very recently Duarte et 
al. (2001) modified the technique developed by Prelle & Singer (1983) and applied 
it to second order ODEs. Their approach was based on the conjecture that if an 
elementary solution exists for the given second order ODE then there exists at least 
one elementary first integral I{t, x, x) whose derivatives are all rational functions of 
t, X and X. For a class of systems these authors (Duarte et al. 2001) have deduced 
first integrals and in some cases for the first time through their procedure. 

In this paper we show that the theory of Duarte et al. (2001) can be extended in 
different directions to isolate even two independent integrals of motion and obtain 
solutions explicitly. In the earlier study it has been shown that the theory can be 
used to derive only one integral. In this work we extend their theory and deduce 
general solution from the first integral. Our examples include those considered in 
their paper and certain important equations discussed in the recent literature whose 
solutions are not known. Our study was motivated by two reasons. Firstly, it is to 
show that apart from finding first integrals one can also deduce general solutions in 
a straightforward and simple manner. Here the method we propose is not confined 
to the PS method alone but can be treated as a general one, that is, suppose one 
has a first integral for a given second order ODE then our method provides the 
general solution in an algorithmic way at least for a class of equations. The reason 
for merging our procedure with PS method rather than any other method is due 
to the following facts. (1) For a given problem if the solution exists it has been 
conjectured that the PS method guarantees to provide first integrals. (2) The PS 
method not only gives the first integrals but also the underlying integrating factors, 
that is, multiplying the equation with these functions we can rewrite the equation 
as a perfect differentiable function which upon integration gives first integrals in 
a separate way. (3) The PS method can be used to solve nonlinear as well as 
linear second order ODEs. As the PS method is based on the equations of motion 
rather than Lagrangian or Hamiltonian the analysis is applicable to deal with both 
Hamiltonian and non-Hamiltonian systems. 

Our second reason is to bring out a novel and straightforward way to construct 
linearizing transformations. Particularly we demonstrate that using our procedure 
one can also explore linearizing transformations in a simple and straightforward 
manner. The latter can be used to transform the given second order nonlinear 
ODEs to a linear one, in particular, to the free particle equation. As we illustrate 
below these transformations can be deduced from the first integral itself which is 
totally a new technique in the current literature. In a nutshell, once a first integral 
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is known then our procedure, at least for a class of problems, not only gives the 
general solutions but also provides linearizing transformations. The ideas proposed 
here can be applied to coupled system of second order ODEs as well as higher order 
ODEs, which will be presented separately. 

The paper is organized as follows. In the following section we briefly describe 
the PrcUc-Singcr method applicable for second order ODEs and indicate certain 
new features in finding the integrals of motion. In §3, we have extended the theory 
in three different directions which indicates the novelty of the approach. The first 
significant application is that the second integral can be deduced straightforwardly 
from the method itself in many cases. The second one is to deduce the general 
solution from the first integral. Finally we propose a method to identify linearizing 
transformations. We emphasize the validity of the theory with several illustrative 
examples arising in different areas of physics in §4. In §5, we demonstrate the 
method of identifying linearizing transformation with three examples including the 
one studied in the recent literature. We present our conclusions in §6. 



2. Prelle-Singer method for second order ODEs 

In this section, we briefly discuss the theory introduced by Duarte et al. (2001) for 
second order ODEs and extend it suitably such that the general solutions can be 
deduced from the modiflcations. Let us consider the second order ODEs of the form 

x=^, P,QGC[t,x,x], (2.1) 

where over dot denotes differentiation with respect to time and P and Q are polyno- 
mials in t, X and x with coefficients in the field of complex numbers. Let us assume 
that the ODE (2.1) admits a first integral I{t,x,x) = C, with C constant on the 
solutions, so that the total differential gives 

dl = Iidt + I^dx + I±dx = 0, (2.2) 

where the subscript denotes partial differentiation with respect to that variable. 
Rewriting equation (2.1) in the form ^dt—dx = and adding a null term S{t, x, x)x 
dt — S{t, x, x)dx to the latter, we obtain that on the solutions the 1-form 

+ Si^ dt - Sdx -dx = 0. (2.3) 

Hence, on the solutions, the 1-forms (2.2) and (2.3) must be proportional. Multi- 
plying (2.3) by the factor R{t,x,x) which acts as the integrating factors for (2.3), 
we have on the solutions that 

dl = R{(t> + Sx)dt - RSdx - Rdx = 0, (2.4) 

where <j) = P/Q. Comparing equations (2.2) with (2.4) we have, on the solutions, 
the relations 

It = R{<i> + xS), 

irj. = —RS, 

h = -R. (2.5) 
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Then the compatibihty conditions, Itx = Ixt, hx = Ixti Ixx = Ixx, between the 
equations (2.5), require that 

D[S] = -,^x+S<Px + S^ (2.6) 
D[R] = -R{S + cj>x), (2.7) 
Rx = RxS + RSx, (2-8) 

where 

D — ^ + X ^ + 

dt dx dx 

Equations (2.6)-(2.8) can be solved in the following way. Substituting the given 
expression of (j) into (2.6) and solving it one can obtain an expression for S. Once 
5* is known then equation (2.7) becomes the determining equation for the function 
R. Solving the latter one can get an explicit form for R. Now the functions R and 
S have to satisfy an extra constraint, that is, equation (2.8). Once a compatible 
solution satisfying all the three equations have been found then the functions R 
and S fix the integral of motion /(t, x, x) by the relation 

I{t, x,x) = j R{^ + xS)dt - j (^RS + ^ j R{(t> + xS)d?j dx 

~ / i fi(0 + xS)dt - J (rS +^ j R{'t> + xS)dt^ dx 



dx. 
(2.9) 



Equation (2.9) can be derived straightforwardly by integrating the equations (2.5). 
Note that for every independent set (S,R), equation (2.9) defines an integral. 

Thus two independent sets, (Si,Ri), i = 1,2, provide us two independent in- 
tegrals of motion through the relation (2.9) which guarantees the integrability of 
equation (2.1). We noticed that since we are solving equations (2.6)-(2.7) first and 
check the compatibility of this solution with equation (2.8), one often meets the sit- 
uation that all the solutions which satisfy equations (2.6)-(2.7) need not satisfy the 
constraint (2.8) since equations (2.6)-(2.8) constitute an overdetermined system for 
the unknowns R and S. In fact, for a class of problems one often gets a set (5*1, Ri) 
which satisfies all the three equations (2.6)-(2.8) and another set {S2,R2) which 
satisfies only the first two equations and not the third, namely, (2.8). In this situ- 
ation, we find the interesting fact that one can use the first integral, derived from 
the set (Si,Ri), and deduce the second compatible solution (5*2. -R2). For example, 
let the set (5*2, i?2) be a solution of equations (2.6)-(2.7) and not of the constraint 
equation (2.8). After examining several examples we find that one can make the set 
(S'2, i?2) compatible by modifying the form of R2 as 

R2=F{t,x,x)R2, (2.10) 

where R2 satisfies equation (2.7), so that we have 

{Ft + xFx + ct>Fx)R2 + FD[R2] = -FR2{S2 + <l>x)- (2.11) 

Further, if is a constant of motion (or a function of it), then the first term on the 
left hand side vanishes and one gets the same equation (2.7) for R2 provided F is 
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non-zero. In other words, whenever i*" is a constant of motion or a function of it then 
the solution of equation (2.7) may provide only a factor of the complete solution 
R2 without the factor F in equations (2.10). This general form of i?2 with S2 can 
now form a complete solution to the equations (2.6)-(2.8). In a nutshell we describe 
the procedure as follows. First we determine S and R from equations (2.6)-(2.7). 
If the set {S, R) satisfies equation (2.8) then we take it as a compatible solution 
and proceed to construct the associated integral of motion. On the other hand if it 
does not satisfy (2.8) then we assume the modified form R2 = F(7i)i?2, where I\ 
is the first integral which has already been derived through a compatible solution, 
and find the explicit form of F(^Ii) from equation (2.8), which in turn fixes the 
compatible solution (^2, R2)- This set (52, R2) can be utilized to derive the second 
integral. 



3. Generalization 

(a) Identifying a second integral of motion 

Duarte et al. (2001) have considered certain physically important systems and 
constructed first integrals. Further, they mentioned that applying the original PS 
algorithm to these first integrals (by treating them as first order ODEs) one can 
deduce the general solution. An interesting observation we make here is that there is 
no need to invoke the original PS procedure in order to deduce the general solution. 
In fact, as we show below, the general solution can be derived in a self-contained 
way. As the motivation of Duarte et al. (2001) was to construct only the first integral 
they reported only one set of solution (S, R) for the equations (2.6)-(2.8). However, 
we have observed that an additional independent set of solution, namely, (S'2,i?2), 
of equations (2.6)-(2.8) may lead to another integral of motion, I2, and if the latter 
is an independent function of Ii then one can write down the general solution for the 
given problem from these two integrals alone straightforwardly. Now the question is 
whether one will be able to find a second pair of solution for the system (2.6)-(2.8) 
and construct I2 through the relation (2.9). After investigating several examples 
we observed the following, (i) For a class of equations, including that of harmonic 
oscillator, equation coming from general relativity and the generalized modified 
Emden equation with constant external forcing one can construct a second pair of 
solution (S'2,i?2) and deduce I2 through the relation (2.9) straightforwardly. We 
call this class as Type I. (ii) For another class of equations we can find explicitly 
(S'2,i?2) from (2.6)-(2.8) but one is unable to integrate equation (2.9) exactly and 
obtain the second integral I2 explicitly. We call this class as Type II. The examples 
included in this category are Helmhotz oscillator and Duffing oscillator. For this 
class of equations we identify an alternate way to derive the second integration 
constant, (iii) Besides the above there exists another category in which the systems 
do not even admit a second pair (^2, R2) of solution in simple rational forms for the 
equations (2.6)-(2.8) and we call this category as Type III. An example under this 
class is the Duffing - van der Pol oscillator which is one of the prototype examples 
for the study of nonlinear dynamics in many branches of science. For this class of 
equations also we identify an alternate way to obtain the second integral. 
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(b) Method of deriving general solution 

To overcome the difficulties in constructing the second constant in Types II and 
III we propose the following procedure. As our aim is to derive the general solution 

for the given problem, wc split the functional form of the first integral / into two 
terms such that one involves all the variables {t, x, x) while the other excludes x, 
that is, 

I = Fiit,x,x) + F2{t,x). (3.1) 

Now let us split the function Fi further in terms of two functions such that Fi 
itself is a function of the product of the two functions, say, a perfect differentiable 
function ^Gi{t,x) and another function G2{t,x,x), that is, 

'^G,{t,x)]+F2{G,{t,x)). (3.2) 



G2{t, X, x) dt 



We note that while rewriting equation (3.1) in the form (3.2), we demand the 

fimction F2{t.x) in (3.1) automatically to be a fimction of Gi{t.x). The reason 
for making such a specific decomposition is that in this case one can rewrite equa- 
tion (3.2) as a simple first order ODE for the variable G\ (see equation (3.4) below). 
Actually, wc realized originally such a possibility for the case of the integrable force- 
free Duffing- van der Pol oscillator equation (Chandrasekar et al. 2004), which has 
now been generalized in the present case. Now identifying the function Gi as the 
new dependent variable and the integral of G2 over time as the new independent 
variable, that is. 



w • 



Gi{t,x), z = I G2{t',x,x)dt', (3.3) 

J o 



one indeed obtains an explicit transformation to remove the time dependent part in 
the first integral (2.9). We note here that the integration on the right hand side of 
(3.3) leading to z can be performed provided the function G2 is an exact derivative 
of t, that is, G2 = ■^z{t, x) = xz^ + Zt, so that z turns out to be a function t and 
X alone. In terms of the new variables, equation (3.2) can be modified to the form 

I = F^(^^'^+F2{w). (3.4) 

In other words, 

Fi(^)=I-F2{w). (3.5) 



dz _ 

Now rewriting equation (3.4) one obtains a separable equation 

f = /(»), (3.6) 

which can lead to the solution after an integration. Now rewriting the solution in 
terms of the original variables one obtains a general solution for (2.1). 



Article submitted to Royal Society 



Integrability and Linearization 



7 



(c) Method of deriving linearizing transformations 

Finally, one may note the following interesting point in the above analysis. 
Suppose F2{iu) is zero in equation (3.4). Then one obtains the simple equation 

where / is a constant. In other words we have 

d^w 

^ = 0, (3.8) 

which is nothing but the free particle equation. In this case the new variables z 
and w helps us to transform the given second order nonlinear ODE into a linear 
second order ODE which in turn leads to the solution by trivial integration. The 
new variables z and w turn out to be the linearizing transformations. We discuss 
in detail about this possibility separately in §5. 



4. Applications 

In this section, wc demonstrate the theory discussed in the previous section with 
suitable examples. In particular, we consider several interesting examples, including 
those considered in Duarte et al. (2001) and derive general solutions and establish 
complete integrability of these dynamical systems. We split our analysis into three 
categories. In the first category we consider examples in which the I^s, i = 1,2, 
can be derived straightforwardly from the relation (2.9). In the second and third 
categories we follow our own procedure detailed in §3 b and §3 c and deduce the 
second constant. We note that one can apply our procedure to a wide range of 
systems of second order of type (2.1) but for illustration purpose we consider only 
a few examples in the following. 

(a) Type - I Systems 

As mentioned earlier for certain equations one can also get the second pair of 
solution (52, R2), in an algorithmic way, from the determining equations (2.6)-(2.8) 
and construct I2 through the relation (2.9). We observe that the Examples 1 and 
2 discussed in Duarte et al. (2001) can be solved in this way and so we consider 
these two examples first and then a nontrivial example in the following. 

Example 1 : An exact solution in general relativity 

Duarte et al. (2001) considered the following equation, which was originally 
derived by Buchdahl (1967) in the theory of general relativity, 

XX = 2,x^ + (4.1) 

and deduced the first integral / through their procedure. In the following, we briefly 
discuss their results and then illustrate our ideas. Substituting (j) = ^ + | into 
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(2.6)-(2.8) we get 



5. + xS. + m^s, = % + ^5 + S\ (4.2) 
tx x^ tx 



R, + xRx + --^±Ar, = -RS - ^-^R, (4.3) 
tx tx 

Rx = SRi + RSi. (4.4) 

As mentioned in Sec. 2 lot us first solve equation (4.2) and obtain an explicit form 
for the function S. To do so Duarte et al. (2001) consider an ansatz for S of the 
form 

g ^ a{t,x) + b{t,x)x 
c{t, x) + d{t, x)x ' 

where a, 6, c and d are arbitrary functions of t and x. A rational form for S can be 
justified, since from (2.5) one may note that 5 = So we consider only rational 
forms in x for 5 for all the examples which we consider in this paper. One may 
note that in certain examples, including the present one and examples 3 and 5 
(given below), this form degenerates into a polynomial form in x. However, for 
other examples like the examples 2 and 4 below, one requires a rational form such 
as (4.5). To be general, we carry out an analysis with the form (4.5). 

Substituting (4.5) into (4.2) and equating the coefficients of different powers of 
X to zero we get a set of partial differential equations for the variables a, b, c and 
d. Solving them one obtains 

5i=--, S2 = --. (4.6) 

X X 

We note that Duarte et al. (2001) have reported only the expression 5*1, as the 
solution for equation (4.2). However, we find ^2 also forms a solution for (4.2) and 
helps to deduce the general solution. Substituting these forms Si and ^2 into (4.3) 
and solving the latter one can get an explicit form for the function R. Let us first 
consider Si. Substituting the latter into (4.3) we get the following equation for R: 

+ iR^ + ^iB^R, = _ ^^R. (4.7) 

In order to solve (4.7) again one has to make an ansatz. We assume that the 
following form for R, that is, 

R = A{t,x) + B{t,x)x, (4.8) 

where A and B are arbitrary functions of (t, x). Since R = —I± (vide equation (2.5)) 
the form of R may be a polynomial or rational in i. So depending upon the problem 
one has to choose an appropriate ansatz. To begin with one can consider simple 
polynomial (in x) for R and if it fails then one can go for rational forms. Let us start 
with equation (4.8). Now substituting (4.8) into (4.7) and equating the coefficients 
of different powers of x to zero and solving the resultant equations one can obtain 
finally Ri = j^- Now the solution Si = and Ri = has to satisfy the 
equation (4.4) in order to be a compatible solution, which is indeed true. Once R 
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and S have been found the first integral / can be fixed easily using the expression 
(2.9) as 

A = |r. (4 0) 

One can easily check that Ji is constant on the solutions, that is, ^ = 0. This 
integral has been deduced in Duarte et al. (2001). However, the second expression, 
^2 has been ignored by the authors since the corresponding i?2 coming out of (4.3) 
does not form a compatible solution, that is, it does not satisfy (4.4). But in the 
following we show how one can make it compatible and use it effectively to deduce 
the second integration constant. 

Now substituting the expression ^2 = — | into (4.3) and solving it in the same 
way as outlined in the previous paragraph we obtain the following form for R, 
namely, 

i?2 = ^. (4.10) 

However, this set (S'2,i?2) does not satisfy the extra constraint (4.4). In fact, not 
all forms of R coming out from (2.7) satisfy (2.8) in general. As we explained in 
§3, the form of R^ given in (4.10) may not be the 'complete form' but might be a 
factor of the complete form. To recover the complete form of R one may assume 
that 

R = F{h)R, (4.11) 

where F{Ii) is a function of the first integral /i, and determine the form of F{Ii) 
explicitly. For this purpose we proceed as follows. Substituting 

^2 = F(/i)i?2 = ^^(^i) (4-12) 
into equation(4.4), we obtain an equation for F, that is, 

/iF' + 2i^ = 0, (4.13) 

where prime denotes differentiation with respect to Ii. Upon integrating (4.13) we 
get (after putting the constant of integration to zero) 

(4.14) 



It 



which fixes the form of i?2 as 



Now one can easily check that this set S'2 = — f and i?2 = |§ is a compatible 
solution for the equations (4.2)-(4.4). Substituting ^2 and R2 into (2.9) we can get 
an explicit form for I2, namely, 

l2 = t{t+^). (4.16) 

X 
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From the integrals Ii and I2 one can deduce the general solution directly (without 
performing any further integration) for the problem in the form 



Of course the same result can be obtained solving equation (4.9) from the first 
integral. However, the point we want to emphasize here is that an independent 
second integral of motion can be deduced to find the solution without any further 
integration, which can be used profitably when the expression for 7i cannot be 
solved straightforwardly. 

Example 2 : Simple harmonic oscillator 

To illustrate the above procedure also works for linear ODEs, in the following 
we consider the simple harmonic oscillator and derive the general solution. As the 
procedure of deriving the first integral has been discussed in detail in Duarte et 
al. (2001) we omit the details and provide only the essential expressions in the 
following. 

The equation of motion for the simple harmonic oscillator is 

x = -x (4.18) 

so that the equations (2.6)-(2.8) become 

St + xSx-xS± = l + S'^, (4.19) 
Rt + xRx - xRi = -RS, (4.20) 
Rx - SR± - RSi = 0. (4.21) 

As shown in Duarte et al. (2001) a simple solution for the equations (4.19)-(4.21) 
can be constructed of the form 

Si = ^, Ri=x (4.22) 
x 

which in turn gives the first integral 

Ji=i;2+x2, (4.23) 
through the relation (2.9). However, one can easily check that 

82 = --, R2 = x (4.24) 

X 

is also a solution for the set (4.19) and (4.20) (which has not been reported earlier) 
but does not satisiy the extra constraint (4.21). So, as before, let us seek an R2 of 
the form 

R2^F{h)R2=F{h)x, (4.25) 

where F{Ii) is a function of Ii. Substituting (4.25) into equation (4.21) and inte- 
grating the resultant equation, we get F = j^. Thus R2 becomes 

R2 = Y = (4.26) 
h x^ + x^ 
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Now one can check that (5*2,1^2) satisfies all three equations (4.19)-(4.21) and fur- 
nishes the second integral through the relation (2.9) of the form 

h = -t- I .^^ ^ dx - I ( . " / -2 ? 2 ^^V^' 

J x'^ + x'^ J \x^ + x"^ dx J x^ + x'^ J 

= -t-tan-^-. (4.27) 
x 

Using (4.23) and (4.27) we can write down the general solution for the simple 
harmonic oscillator directly in the form 

x = y/hcos{t + h). (4.28) 

In a similar way one can deduce general solution for a class of physically important 
systems. 

One may note that in the above two examples, I2 can also be obtained trivially 
by simply integrating the expressions (4.9) and (4.23) without using the extended 
procedure. We stress that for certain equations one is not able to integrate and 
obtain the general solution in this simple way and has to follow the above said 
procedure in order to obtain the second integral. In the following we discuss one 
such example for which to our knowledge explicit solution was not known earlier. 



Example 3 : Modified Emden type equation with linear term 

It is known that the generalized Emdcn type equation with linear and constant 
external forcing is also linearizable since it admits an eight point Lie symmetry 
group (Mahomed & Leach 1989; Pandey et al. 2004). In the following we explore 
its general solution through the extended PS algorithm. Let us first consider the 
equation of the form 

x + kxx + —x^ + \ix = 0, (4.29) 

where k and Ai are arbitrary parameters. To explore the general solution for the 
equation (4.29) we again use the PS method. In this case we have the following 
determining equations for the functions R and 5, 

St + xS^ - {kxx + —x^ + \ix)Si =kx + —x^ + Ai - Skx + 5^ (4.30) 
y o 

k'^ 

Rt + xRx - {kxx + —x^ + \ix)R± = -R{S - kx), (4.31) 
9 

R^ - SRi - RSi = 0. (4.32) 

As before, let us seek an ansatz for S of the form (4.5) to the first equation in 
(4.30)-(4.32). Substituting the ansatz (4.5) into (4.30) and equating the coefficients 



Article submitted to Royal Society 



12 



Chandrasekar, Senthilvelan and Lakshmanan 



of different powers of x to zero we get 



dbx — bdx — 


kd^ 


= 0, 


dbt — bdt + cbx — bcx + axd — adx — 2kcd — {-^x^ + Ai)rf^ + kbdx 

O 


-b^ 


= 0, 


cbt — bct + dat — adt + ca^ — acx — kc^ — 2{—x^ + \i)cd + 2kadx — 


2ab 


= 0, 


P k^ 
cat — o-ct — {~frx^ + Mx){bc — ad) — (-^a;^ + Ai)c^ + kacx ■ 




= 0, 






(4.33) 



where subscripts denote partial derivative with respect to that variable. Solving 
equation (4.33) we can obtain two specific solutions, 

-X + kx + 3y/^ kx 

^^ = ^^' = 3 kx + 3V^,- ^"-^'^ 

Putting the forms of 5i and 5*2 into (4.31) and solving it one can obtain the re- 
spective forms of R. To do so let us first consider Si. Substituting the latter into 
(4.31) we get the following equation for R: 

Rt + xRx - {kxx + —x^ + \ix)Rx = ( \-kx\R. (4.35) 

Again to solve equation (4.35) we make an ansatz of the form 

R = A{t,x)+B{t,x)x 

C{t,x) + D{t,x)x + E{t,x)x^' ^ ' 

Now substituting (4.36) into (4.35) and equating the coefficients of different powers 
of X to zero and solving the resultant equations we arrive at 

where Co = 18\/— Ai. One can easily check that Si and Ri satisfies equation (4.32) 
and as a consequence one obtains the first integral 

I,=e-^^-('^^±^±^^). (4.38) 
\3x + kx'' — 3v— Aix J 

We note that unlike the other two examples, equation (4.38) cannot be inte- 
grated straightforwardly to provide the second integral (though one can in fact 
explicitly solve the resultant Riccati equation after some effort). We follow the pro- 
cedure adopted in the previous two examples and construct 12- Now substituting 
the expression 5*2 into (4.31) and solving it in the same way as outlined above we 
obtain the following form for R, that is, 

k{3x + kx^ - 3^/^x)^ ^ ' 
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However, this set (S'2,i?2) does not satisfy the extra constraint (4.32) and so to 
deduce the correct form of R2 we assume that 



6 -FfMP _^ -F(Ji)(fca: + 3V^)e-3v^* 

R, _ F{I,)R2 - Co fc(3^ + fc^2_3^^)2 • (4-40) 

Substituting (4.40) into equation (4.32) we obtain F = j^, which fixes the form of 
R as 



Now one can easily check that this set (^2, -R2) is a compatible solution for the set 
(4.30)-(4.32) which in turn provides I2 through the relation (2.9), 

2 ,/3A7t / 9Ai + 3kx + k^x^ 



y-ix + kx"^ + ■i^/^x)' ^^'^^^ 

Using the explicit form of the first integrals I\ and I2 , the solution can be deduced 
directly as 



3V^(/ie2^=^* - 1) 
khhe^^t + k{l + 7ie2v^t) J ' 



(4.43) 



To our knowledge, the above explicit solution, (4.43), of the equation (4.29) is given 
for the first time. It has several interesting consequences for nonlinear dynamics, 
which will be discussed separately. 



(6) Type - // Systems 

In the previous category we considered examples which straightforwardly give 
the integrals Ii and I2 through the relation (2.9). In the present category we show 
that there are situations in which an explicit form of I2 is difficult to obtain through 
the relation (2.9) even though one has a compatible solution for (2.6)-(2.8). So one 
has to go for an alternate way in order to obtain the general solution for the given 
problem. For this purpose, we make use of the method proposed in §3 6. In the 
following we give examples where such a possibility occurs and how to over come 
this situation. 



Example 4 ■' Helmholtz oscillator 

Recently Almendral & Sanjuan (2003) studied the invariance and integrability 
properties of the Helmholtz oscillator with friction, 

X + C1X + C2X - f3x'^ = 0, (4.44) 

where Ci, C2 and /3 are arbitrary parameters, which is a simple nonlinear oscilla- 
tor having a quadratic nonlinearity. Using the Lie theory for differential equations 

6c^ 

Almendral & Sanjuan (2003) found a parameteric choice C2 = for which the 
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system is integrable and derived the general solution for this parametric value. In 
the following we solve this problem through the extended PS method. 

Substituting </> = — {cix + 02.x — fix^) into equations (2.6)-(2.8) we obtain 

St + xS^ - {cix + cix - l3x^)Si = C2 - 2(3x - ciS + S"^, (4.45) 
Rt + xR^ - {cix + C2X - (3x^)Ri = -R{S - a), (4.46) 

R^ = SRi + RSi. (4.47) 



Making the same form of an ansatz, vide equations (4.5) and (4.8), we find nontrivial 

25 



6c 

solution exists for (4.45)-(4.46) only for the parametric restriction C2 = The 



respective solutions are 



S^^L^^^I^^^ i?i = -(i+^a.)et-*, (4.48) 

X + ^x 5 

82=^-^—^ i?2 = -i;e^^*. (4.49) 

X 

Now one can easily check that {Si, Ri) satisfies the third equation (4.47) and as a 
consequence leads to the first integral of the form 

e^^^fx^ ^ 2cixx ^ 2c\x'^ I3x^ 



However, the second set (S'2,-R2) does not satisfy the extra constraint (4.47) and 
so we take 

R2 = F{I)R2 = -F{I)xe'''\ (4.51) 
which in turn gives F = CoI~i , where Co is an integration constant, so that 

\lfJ (f + ^ + ^-/3f)t 

One can check that (5*2, -R2) satisfy equations (4.45)- (4.47) and so one can proceed 
to deduce the second integration constant through the relation (2.9). However, upon 
substituting {S2,R2) into (2.9) we arrive at 



■dx. (4.53) 



. . 6Ci X r, 2 

cix + px 

V 2 5 """25 Z-* 3 / 

It is very difficult to evaluate the integral and so is one is not able to obtain an 
explicit form of I2 for this problem through this way. A similar form of I2 has been 
also derived by Bluman & Anco (2002) and Jones et al. (1993) for the Duffing 
oscillator problem (that is cubic nonlinearity in equation (4.44)). 

Unlike the other examples discussed in Type I the present example possesses 
difficulties in evaluating the second integration constant. In fact, for a class of equa- 
tions one faces such complicated integrals. To overcome this one has to look for an 
alternate way such that the second constant can be deduced in a straightforward 
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and simple way. We tackle this situation in the following way. As we have seen, in 
most of the problems, we are able to deduce the first integral, that is, 7i, straightfor- 
wardly and the first integral often admits explicit time dependent terms. A useful 
way to overcome this is to remove the explicit time dependent terms by trans- 
forming the resultant differential equation into an autonomous form and integrate 
the latter and obtain the solution. In order to do this one needs a transformation 
and the latter can often be constructed through ad- hoc way. However, as we have 
shown in the theory in §3 6 , one can deduce the required transformation coordinates 
in a simple way from the first integral itself and the problem can be solved in a 
systematic way. 

Rewriting the first integral Ii given by equation (4.50), in the form (3.1), we 

get 

^ " 2('^"^~5~J ^ 3 ^ ^ ^ 

Now splitting the first term in equation (4.54) further in the form (3.2), 

/i=e^(|(^xe^))'-f(xei-T, (4-55) 

and identifying the dependent and independent variables from (4.55) and the rela- 
tions (3.3), we obtain the transformation 

1 2cit 5 _ £l£ 

w = —=xe 5 , z = e 5 . (4.56) 

V2 ci ^ ^ 

One can easily check that equation (4.44) can be transformed to an autonomous 
form with the help of the transformation (4.56). Wc note that the transformation 
(4.56) exactly coincides with the earlier one which has been constructed via Lie 
symmetry analysis in Almendral & Sanjuan (2003). 

Using the transformation (4.56) the first integral (4.54) can be rewritten in the 
form 

/ = w'^ - ^w^ (4.57) 

which in turn loads to the solution by an integration. On the other hand the trans- 
formation changes the equation of motion (4.44) to 

w" = Pw^, (4.58) 

where /3 = 2^/2(3, which upon integration gives (4.57). Prom equation (4.57), we 
obtain 

= Aw^ - ga, (4.59) 

where z = 2^/|i and gs — The solution of this differential equation can 

be represented in terms of Weierstrass function g{z; 0, 53) (Gradshteyn & Ryzhik 
1980; Almendral & Sanjuan 2003). 
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(c) Type - III Systems 

In the previous two categories we met the situation in which we are able to con- 
struct a pair of solutions ^2) for the equations (2.6) from which Ri and R2 have 
been deduced. However, there arc situations in which one is able to construct only 
one set of solution ^i) and its corresponding first integral only and the second 
pair of solution {R2, S2) can not be obtained by simple rational form of ansatz. In 
this situation one can utilize our procedure and deduce the general solution for the 
given problem. In the following we illustrate this with a couple of examples. 



Example 5 : Force free Duffing-van der Pol oscillator 

One of the well-studied but still challenging equations in nonlinear dynamics is 
the Duffing-van der Pol oscillator equation. Its autonomous version (force-free) is 

x + {a + px'^)x -^x + x^ = 0, (4.60) 

where over dot denotes differentiation with respect to time and a, (3 and 7 are 
arbitrary parameters. Equation (4.60) arises in a model describing the propagation 
of voltage pulses along a neuronal axon and has received a lot of attention recently 
by many authors. A vast amount of literature exists on this equation, for details see 
for example Lakshmanan & Rajasekar (2003) and references therein. In this case 
we have 

St + xSx - ((a + (3x'^)x -^x + x^)Si = {2(5xx - 7 Sx^) 

-{a + l3x^)S + S'^, (4.61) 
Rt + xRx - {{a + j3x^)x - jx + x^)R± = {a + (3x^ - S)R, (4.62) 

Rx = SRi+RSi. (4.63) 

To solve equation (4.61)-(4.63) we seek an ansatz for S and R of the form 

_ a(t, X) + b(t, X)X „ ^ T,/ ^. , . „ 

'=W-^;kw±x^ i? = ^(M) + 5(t,.).. (4.64) 



Upon solving the equations (4.61)-(4.63) with the above ansatz we find nontrivial 



solution exists only for the choice a = 4 , 7 = — and the corresponding forms of 



S and R reads 

S=^+f3x'^, R = eT. (4.65) 

For this set one can construct an invariant through the expression (2.9) which turns 
out to be (Senthilvelan & Lakshmanan 1995) 

Id 3 
x + -x+^x^ = Ie~^K (4.66) 

To obtain a second pair of solutions for the equations (4.61)-(4.63) one may seek 
more general rational form of S and R by including higher polynomials in x. How- 
ever, they all lead to only functionally dependent integrals. As it is not possible to 
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seek the second pair of solution by simple ansatz one has to see an alternate way as 
indicated in §3 6 . We can deduce the required transformation coordinates from the 
first integral and transform the latter to an autonomous equation and integrate it. 

Using our algorithm given in §3 6 one can deduce the transformation coordinates 
from the first integral itself which turns out to be (Chandrasekar et al. 2004) 



w = -xeT^ \ z = e ^ , (4.67) 

where w and z are new dependent and independent variables respectively. Sub- 

^ -7 = -4 



stituting (4.67) into (4.60) with the parametric restriction a = 4,7 = —4, we 



get 



02 

w" - ^v?w' = 0, (4.68) 



where prime denotes differentiation with respect to z. Equation (4.68) can be inte- 
grated trivially to yield 

w' - y = /, (4.69) 

where I is the integration constant. Equivalently, the transformation (4.67) reduces 
(4.60) to this form. Solving (4.69), we obtain (Gradshteyn & Ryzhik 1980) 



a 

'-''=31 



1 / (w + a)2 \ /- / wV^ 
o ^ + \/3arctan 

2 yw'' — aw + a'' J \2a — w 



(4.70) 



where <i = \J jj^ a,nd zo is the second integration constant. Rewriting w and z in 

terms of old variables one can get the explicit solution for the equation (4.60). 

In the above, we have shown that the systems (4.44) and (4.60) are integrable 
for certain specific parametric restrictions only. One may also assiimc that the 
functions S and R involve higher degree rational functions in x and repeat the 
analysis. However, such an analysis does not provide any new integrable choice. In 
fact, the present results coincide exactly with the results obtained through other 
methods, namely, Painleve analysis, Lie symmetry analysis and direct methods 
(Senthilvelan & Lakshmanan 1995; Almendral & Sanjuan 2003; Lakshmanan & 
Rajasekar 2003). 



5. Linearizable equations 

In the previous section we discussed the complete integrability of nonlinear dynam- 
ical systems by constructing sufficient number of integrals of motion and obtaining 
the general solutions explicitly. Another way of solving nonlinear ODEs is to trans- 
form them to linear ODEs, in particular to a free particle equation and explore their 
underlying solutions. Eventhough this is one of the classical problems in the theory 
of ODEs, recently considerable progress has been made (Mahomed & Leach 1989a; 
Steeb 1993; Olver 1995; Harrison 2002). In this direction it has been shown that a 
necessary condition for a second order ODE to be linearizable is that it should be 
of the form (Mahomed & Leach 1989a) 

q = D{t, q) + C{t, q)q + B{t, q)q^ + A{t, q)q^, (5.1) 
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where the functions A, B, C, and D are analytic. SufScient condition for the above 
second order equation to be Unearizable is (Mahomed & Leach 1989a), 

3Att + 3CAt - 3DAg + 3ACt + Cgq - 6ADg + BCg - 2BBt - 25(, = 0, 

Btt + &DAt - 3DBq + 3ADt - 2Ctq - 3BDg + 3Dgg + - CBt = 0, (5.2) 

where the suffices refer to partial derivatives. 

For a given second order nonlinear ODE one can easily check whether it can be 
Unearizable or not by using the above necessary and sufficient conditions. However, 
the nontrivial problem is how to deduce systematically the linearizing transforma- 
tions if the given equation is linearizable. As far our knowledge goes Lie symmetries 
are often used to extract the linearizing transformations (Mahomed & Leach 1985). 
As we pointed out in §3 the linearizing transformations can also be deduced from 
the first integral itself, whenever the system is linearizable, in a simple and straight- 
forward way and we stress that our procedure is new to the literature. In fact, we use 
the same procedure discussed in §3 c and deduce the linearizing transformations. 
The only difference is that in the case of linearizing transformations the function 
F2 turns out to be zero in equation (3.2) and as a consequence the latter becomes 
^ = / and the transformation coordinates become the linearizing transformations. 
We illustrate the theory with certain new examples in the following. 



Example 1 : General relativity 

To illustrate the underlying ideas let us begin with a simple and physically 
interesting example, namely, the general relativity equation which we discussed as 
Example 1 in §4. We derived the solution (4.17) using the PS method. In the present 
section we linearize the system and derive its solution. Rewriting the first integral 
(4.9) in the form (3.1) 

and identifying (5.3) with (3.2), we get 

Gi = ^, G2 = -2t, F2=0. (5.4) 

With the above choices, equation (3.3) furnishes the transformed variables, 

w = \, z = -t^. (5.5) 

Substituting (5.5) into (4.1), the latter becomes the free particle equation, namely, 
,2 

= 0, whose general solution is w = hz + I2, where h and I2 are integration 
constants. Rewriting w and z in terms of x and t one gets exactly (4.17) which has 
been derived in a different way. 



Example 2 : Modified Emden type equations 

Recently several papers have been devoted to explore the invariance and inte- 
grable properties of the modified Emden type equations (Mahomed & Leach 1985; 
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Duarte et al. 1987), 



X + kxx + —x^ = 0. 



(5.6) 



In fact, it is one of the rare second order nonlinear ODEs which admit eight Lie 
point symmetries and as a consequence is a linearizable one. Recently Pandey et 
al. (2004) have obtained tlic' explicit forms of the Lie point symmetries associated 
with the more general equation 



where k, Ai and A2 are arbitrary parameters. They found that not only the Emden 
equation (5.6), but also its general form, that is, equation (5.7), admits eight Lie 
point symmetries. The authors have also reported that the explicit forms of the 
symmetry generators. However, due to the complicated forms of the symmetry 
generators it is difficult to derive the first integrals and linearizing transformations 
from the symmetries straightforwardly (though in principle this is always possible) . 
Nevertheless, we discussed about the integrability of the case A2 = 0, Ai of 
equation (5.7) as Example 3 in §4 and deduced its general solution. In this section we 
transform the equation into free particle equation and deduce the general solution in 
an independent way. We divide our analysis into 2 cases, namely, (i) Ai ^ 0, A2 = 
and (ii) Ai 7^ 0, A2 7^ and construct linearizing transformations and general 
solutions for both the cases. As the procedure is same as given in the previous 
examples we give only the results. 

Case(i) A2 = 0, Ai ^ ; Modified Emden type equation with linear term 
Restricting A2 = in (5.7) we have 



Since the first integral is already derived, vide equation (4.38), we utilize it here 
to deduce the linearizing transformations. Rewriting the first integral (4.38) in the 
form 



X + kxx + ^x^ + X\x + A2 = 



(5.7) 



X + kxx + —x^ + Aia; = 0. 



(5.8) 




(5.9) 



and identifying (5.9) with (3.2), we get 




(5.10) 



With the above functions (3.3) furnishes 




(5.11) 
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which is nothing but the hnearizing transformation. One may note that in this 
case also while rewriting the first integral / (equation (4.38)) in the form (3.1), the 
function F2 disappears, and as a consequence we arrive at (vide equation (3.4)) 




(5.12) 



which in turn gives the free particle equation by differentiation or leads to the 
solution (4.43) by an integration. On the other hand vanishing of the function F2 
in this analysis is precisely the condition for the system to be transformed into the 
free particle equation. 



Case(ii) Ai ^ 0, A2 7^ ; Modified Emden type equation with linear term and 
constant external forcing 

Finally, we consider the general case, that is, 

X + fcxi; + ya;^ + Aix + A2 = 0. (5.13) 

To explore the first integrals associated with the system (5.13) let us seek the PS 
algorithm again. The determining equations for the functions R and S turn out to 
be 

St + xSx - {kxx + —x^ + Xix + X2)S± =kx+ —x^ + Ai - Skx + S^, (5.14) 
9 3 

Rt + xR^ - {kxx + — a;^ + X-^x + A2)-Ri = {kx - S)R, (5.15) 

R^ - SRi - RSi = 0. (5.16) 

As before let us seek an ansatz for S to solve the equation (5.14), namely, 

a{t, x) + b{t, x)x 



c{t, x) + d{t, x)x 



(5.17) 



Substituting (5.17) into (5.14) and equating the coefficients of different powers of 
X to zero and solving the resultant equations we arrive at 

kx + 3a kx kx + 3/3 kx 

where + aAi - ^ = and /3 = -«±V-3a^-4A i _ pitting the forms of Si into 
(5.15) we get 

Rt + xR^ - {kxx + + Xix + X2)Ri = ( + kx\ R. (5.19) 

Again to solve this equation we make an ansatz 

D = Mt, x)+B{t,x)x 

C{t, x) + D{t, x)x + E{t, x)x^ ■ ^ ' 
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Substituting (5.20) into (5.19) and solving it we obtain the following form of R, 
namely, 



{3kx - S^^^^ikx + 3a) + {kx + 3a)2)2 ' 

where Co is constant and a = V— Sa^ — 4^^. We find that the solution {S\,R{) 
satisfies (5.16). Equations (5.18) and (5.21) fix the first integral of the form 

J ^ / Sfcx - 3i^(fcx + 3a) + {kx + 3a)2 \ ^ 
^ ^v3fci_3(3^(fca; + 3Q,) + (fca; + 3c)2y'' 

where Cq = 9fcd. Rewriting the first integral (5.22) in the form (3.1) 

e=^\kix + iaf 



Zkx - ^^^^^{kx + 3a) + {kx + 3a)2 
d ( . —3 3a ± a , i,<x^& ^ 



.rft V'fca; + 3a ' 2(3a2 + Ai)''' ' ^^'^^^ 

and identifying (5.23) with (3.2), we get 

_ / -3 3a±d \ isiHj 

^ ~ \kx^2,a ^ 2(3a2 + Ai)j'' ' ' 

Zkx-Z''^^^{kx + 'ia) + {kx^2,af 3a±a, 
= (L + 3a)2 (5-24) 



so that (3.3) gives 



—3 3a ± d \ 3QTa + 

u) = I 1 p — 2 

\A;:c + 3a 2(3a2 + Ai)y 

^=(^ + 2(3a2+l) )^"^^' ^'-''^ 

which is nothing but the linearizing transformation. Substituting (5.25) into (5.13) 
we get the free particle equation 

,5.e, 

whose general can be written as w = Iiz + l2- Rewriting w and z in terms of the 
original variable x and t one obtains 



. = + ^ f (3^^ + A0(l-ir^) X ^ 

k k \sa{l - 7ie±«*) - 2(3a2 + Ai)72e^^* ± d(l + 7ie±«*) / 

On the other hand the general solution can also be derived by extending the PS 
method itself. To do so one has to consider the function 82- Now substituting the 
expression ^2 into (5.15) and solving it in the same way as outlined in the previous 
paragraphs we obtain the following form for R, that is, 

, C.(^x + 3ffle'^. 

(3fca; - 3^^^(fca; + 3a) + (fca; + 3a)2)2 
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However, this set, (5'2,i?2), does not satisfy the extra constraint (5.16) and to 
recover the full form of the integrating factor we assume that 

R2 = F{h)R2. (5.29) 

Substituting (5.29) into equation (5.16) we obtain an equation for F, that is, IiF' + 
2F = 0, where prime denotes differentiation with respect to /i. Upon integrating 
this equation we obtain F = l/I^, which fixes the form of R as 

Co(fa + We^- 

{3kx-3^2^{kx + Sa) + {kx + 3a)^)^ 

Now one can easily check that this set S2 and R2 is a compatible solution for the 
equations (5.14)-(5.16). Substituting 52 and .R2 into (2.9), we can obtain an explicit 
form for the second integral I2, that is, 

j__f 2a{3kx-3akx + k'^x'^ + 9a'^ + 9Xi)e^^* \ 
V (3a ± a){3kx - 3^^^(fcx + 3a) + {kx + 3a)^))' 

Rewriting equation (5.22) for x and substituting it into (5.31) we get the same 
expression (5.27) as the general solution. 

Example 3 : Generalized modified Emden type equation 

Recently, Pandey et al. (2004) have considered the following Lienard equation 

x + f{x)x + g{x)=(i, (5.32) 

where / and g are arbitrary function of their arguments, and classified system- 
atically all polynomial forms of / and g which admit eight Lie point symmetry 
generators with their explicit forms. They found that the most general nonlinear 
ODE which is linear in x whose coefficients are functions of the dependent variable 
alone should be of the form 

x + {kix + k2)x + ^x^ + h^x^ + Xix + A2 = 0, (5.33) 

where fcj's and A^'s, i=l,2, arc arbitrary parameters, which is consistent with the 
criteria (5.1) and (5.2) given by Mahomed & Leach (1989a). Interestingly, equa- 
tion (5.33) and all its sub cases posses sl{3, R) symmetry algebra. For example, we 
discussed the intcgrability and linearization of equation (5.33) with ^2 = in the 
previous example. As the linearizing transformations and the general solution of 
equation (5.33) are yet to be reported we include this equation as an example in 
the present work. As in the previous case we divide our analysis into three cases, 
(i) Al = 0, A2 = : Modified Emden type equation with quadratic and cubic 
nonlinearity 

^+(fciar + fc2)i+^x3 + ^a;2 =0. (5.34) 
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(ii) Ai 7^ 0, A2 = : Modified Emden type equation with quadratic and Unear terms 

X + {kix + k2)x + -^x^ + -^x'^ + Xix = 0. (5.35) 

(iii) Ai 7^ 0, A2 7^ : The full generalized modified Emden type equation (5.33). 
We have derived the integrating factors, integrals of motion, linearizing trans- 
formations and the general solutions for all the cases. As the calculations are the 
similar to the one discussed in the previous case, we present the results in tabular 
form (Table I), where the results for the most general case (5.33) has been given 
from which the results for the limiting cases (5.34) and (5.35) can be deduced. 

6. Conclusion 

In this paper we have discussed the method of finding general solutions associated 
with second order nonlinear ODEs through a modified PS method. The method 
can be considered as a direct one, complimenting the well known method of Lie 
symmetries. In particular, we have extended the theory of Duarte et al. (2001) such 
that one can recover new integrating factors and their associated integrals of motion. 
These integrals of motion can be utilized to construct the general solution. In the 
situation when one is not able to recover the second integral of motion we introduced 
another approach to derive the second integration constant. Interestingly, we showed 
that in this case it can be derived from the first integral itself in a simple and 
elegant way. Apart from the above we introduced a technique which can be utilized 
to derive linearizing transformation from the first integral. We illustrated the theory 
with several new examples and explored their underlying solutions. 

In this paper we concentrated our studies only on single second order ODEs. 
In principle the method can also be extended to third order ODEs and systems of 
second order ODEs. The results will be published elsewhere. 
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Table 1. Integrating factors, integrals of motion, linearizing transformations and the 
general solution of equation (5.33) 

Null forms and Integrating factors 

_ A;ia; + 3a kix „ Co(A;ia; + 3a)e'F"* 

~ iti — 



3 A;ia; + 3a (3/fcia; - ^(3feia; + 9a) + (fcia; + 3a)2)2 

kix + 3P kix „ Co{kix + 3l3)e^* 



kix + -i(i' 3fcja;_ m(3jfcia; + 9a) + (ifcia; + 3Q)2 



— k'zQ^ + f-iAi — — ^ =0, a = y^— 3a^ + 2ak2 + fc| — 4Ai, /3 = 3q; — fc2 

„ _ —a + k'z ±a 

P o 5 



First Integrals 

^e,t( ^kix - ^(3fcix + 9a) + (fcia; + 3q)" \ . 

71=6^ I J— I , Go = 9/cia 

V3fcii - ^{?,kix + 9a) + (fcix + 3a) 



/2 



-2ae^* f ^kix - 3fcia:(a - fc2) + fefx^ + 9a^ - 9afc2 + 9Ai 
j3±a V sfeii- ^(3feix + 9a) + (feia: + 3a)2 



Linearizing Transformations 



-3 /3±d \ i3. 



A;ia; + 3a 2(3a2 - 2afe2 + Ai) 



-3 ^ 13 T a 



A;ia; + 3a 2(3a2 - 2afc2 + Ai) 



Solution 



3a_^J_/^ 6(3a2-2afe2 + Ai)(l-/ie=^**) 



ki ki V^(i _ 7ie±««) ± (/3 ± a)/i/2e^^* ± a(l + 7ie±«*) 
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